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1. INTRODUCTION 

Web transport system (WTS), which is one of the popular systems, is widely used in the production 
and handling of web materials [1J—(B] such as paper production, glass, oled materials, solar cells, fabrics, and 
printing operations [4], [5]. The materials have characteristics such as thin, continuous, elastic, and easy to 
be damaged during transportation and handling, so WTS requires not only high technology in terms of me- 
chanics but also high requirements on control engineering. In fact, the controllers used for WTS are mainly 
proportional-integral (PI) and proportional-integral-derivative (PID) [6]-[8], feedforward control [9], and PID 
with feedforward control [10]. However, this is a linear control method, while WTS is a strongly nonlinear sys- 
tem, affected by external noise, so it does not always give the desired quality, but the often slow response, low 
stability, low accuracy, and challenge to meet the increasingly high requirements of today’s modern production 
lines. Therefore, it is necessary to build a modern controller to overcome the disadvantages. 
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In recent years, various nonlinear control methods such as sliding mode control [11-13], intelligent 
control [14-16], backstepping control (BC) [17J-[19], and disturbance compensator have been applied 
to improve the control quality for the system. These control methods can only improve performance when the 
system parameters are clearly known. However, the WTS system is a time-varying system, and the calculation 
of system parameters is quite complicated and less precise. To further improve the system performance, it is 
necessary to have an adaptive mechanism to track the system parameters’ variation and radial basis function 
(RBF) neural network is an appropriate selection in many applications [21-24]. The paper proposes adaptive 
BC for WTS using RBF neural networks. The controller is designed using the backstepping technique while 
the system’s time-varying parts are approximated by RBF neural network. The RBF network weight updating 
rules are proposed such as the system stability and desired system performance are guaranteed. 


2. WEB TRANSPORT SYSTEM DYNAMIC MODELING 

Let us consider a single-span WTS including an unwind roll and a rewind roll as shown in Figure [I] 
The guide rolls are ignored for simplicity. According to Newton’s law and mass conservation, the nonlinear 
dynamic equations of web transport can be written as (I)-(3) (25): 


tw = CyWy + CoW pty + C3Wr (1) 
Wu = C4Mu + C5tw + C6Wu + c7w2 (2) 
Wr = cg Mr + cotw + Ciowr + cw? (3) 


where c;, (i = 1, .., 11) are the time-dependent parameters defined as: 


Nig 7 ESry, Tr ESr, 1 Tu bfu awpr? 1 
C1 = L 3c —4 3 C. = 3 C. = 3C = 3C = 3c = ;C = 3 
PE Ee N GA ee ee io ee 
Tr bfr awp? 
C9 = ; 10 = C11 = ; 
9 J, 10 J, 11 J, 


In the definition of c;, operating radius ru, ry and moment of inertia J,,, Jr, are calculated as: 


1 
rault) = Ruo — Lair, (t) = Rro + La; Jalt) = Juo + spwm(rs, — Rê); 


20 an 
1 
Jn(8) = Jno + Spwn(r — RA) 


Where tọ is web tension, w,, is unwind roll’s angular velocity, w, is rewind roll’s angle velocity, Mu is unwind 
roll torque, M, is rewind roll torque, ru is unwind roll radius, Ruo is unwind roll’s initial radius, r, is rewind 
roll radius, Ro is rewind roll’s initial radius, J„ is the unwind roll’s total moment of inertia, Juo is unwind 
roll’s initial total moment of inertia, J, is rewind roll’s total moment of inertia, Jyo is rewind roll’s initial total 
moment of inertia, bp, is unwind roll vicious friction coefficient, bs, is rewind roll vicious friction coefficient, 
E is web elasticity, S' is web cross-sectional area, L is web length, a is web thickness, w is web width, and p is 
web density. 
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Figure 1. Single-span WTS 
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3. DESIGN A BACKSTEPPING CONTROL INCORPORATING RADIAL BASIS FUNCTION NEU- 
RAL NETWORK 
At first, we define intermediate variables fu, gu, fr, gr as: 


2, = š = 2; _ 
fu = C5tw + C6Wy + CTW; Ju = C4; fr = Cotw + C10Wr + C11W 75 Ir = Cs 


Then, the dynamical model (I}-(3) becomes (4)-(6}: 


t = cyWy, + Cwrt + C3wr (4) 
Wp = fr + grM, (6) 


The designing process is as: 
— Step 1: defining tension tracking error variables as p}: 


Atw = ty — Ta (7) 
where T4 is reference web tension, then: 
Aty =t— Th = C1Wu + CoW,t + C3wWr — T (8) 


Choose the Lyapunov function: Vj = At. Taking the derivative of V; one can obtain as {o}: 


vi = Aight, = (C1wu + CoWpt + C3wr — Ta)(t z Ta) (9) 
In order to stabilize the subsystem, the condition Vi < 0 must be guaranteed. Thus, we choose {10}: 
Cito, + wt + ti, — Ty = —k(t — Ta) (10) 


where k; is a positive real number, then: 


Vi = —ki(t — Ta)? = —k At? < 0; Yki > 0 (11) 
From condition (10), we deduce the virtual control signal to stabilize the subsystem as (12): 


Fug = — E (cztwr + egy — Ta + ke(t — Ta)) (12) 
— Step 2: defining unwind velocity tracking error variables as (13): 


Awu = Wy — Funa (13) 


Choose the Lyapunov function: V, = $(Aw,)? = $ (Wy — Fpa)”. Taking the derivative of V,,: 


1 
2 
Vu = Aw, Aw, = Awy (Wu — Poss) = Awu(fu + GuMu — Pring) (14) 

In order to obtain V, < 0 the control signal M, is chosen as (15): 
Mau = = (fu =F, + kw, Awu) (15) 

where kw, is a positive real number, then: 

Va = —ku, (Awu)? < 0 (16) 
However fu and gu are difficult to determine precisely. Thus, we will approximate them by fa and ĝu, respec- 


tively. The control signal then is calculated as (I7): 
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Ma = -4 (fa — Fig + ku, Awu) (17) 
The derivative of V, becomes (18): 


Vu = Awu((fu- fu) + fu + (Ju = Gu) Mu + GuMu — Pg) 
= Awu(fu = fu) F Awu (Gu = Gu)Mu = Ran, (Awu)? 
Next, the approximation rules for fu and gu are proposed. Supposing that the function f,, can be computed by 
an ideal RBF as (19): 


(18) 


fu = WE tn (19) 
Where W, is the ideal weight vector of the neural network W, = [W Wža, <- Wë m]? m is the number of 
neural in the network; and h, is a Gaussian function vector hy = [hu1, hu2, ---, hum)’. The Gaussian function 


hui, (i = 1,2, ..., m), with inputs wu, wr, is defined as (20): 


12 


2 
wu ~er llf +ller —c9i 
exp( li = 2il”) 


hui = a p 20) 


eu —e1 jl? tllor cz; 
om, exp( u i Ds r i ) 


Where Cii and cz; are the position of the center of the RBF, bu; is the width of the RBF. The approximation 
function fu are calculated using RBF as (Figure h): 


fa = Wa hu (21) 
Where W,, = [Wu Wu2; ---, Wum]? is weight vector. We have to establish the weight updating rule so that the 
system is stable. In addition, we also have to find the approximation rule for gu. It is noted that gu = —1/Ju, 


then there exists upper bound gu < gum < 0. 


Figure 2. Schematic diagram of RBF neural network to approximate fa and T 


Supposing that gu can be approximated by ĝu. We define Wu = Wa — Wa; Ju = Ju — ĝu. Choose 
the Lyapunov candidate function Vz, as (22): 
Vou = Va + WET Wa + img? = 4 (Awu)? + SWI Wa + img? (22) 
where I", is a positive definite diagonal matrix and 7, is a positive real number, then: 


Vou = Va + WETS Wu + Nuĝuĝu . 

= Awul fu = fu) F Awu (gu m ĝu)Mu ku, (Awu)? T wir W, T Nu ĝu ğu 

= Aw WE hy + AwuĝuMu — ku, (Awu)? + WITZ Wu + mududu (23) 
= —kw, (Awu)? + WF (Awuhu +13 Wu) + ĝu ( Awu Mu + Nuĝu) 

= —k,,, (Awu)? + WF (Awuha — rzw.) + Ju (Awu Mu — Nuĝu 
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In order to guarantee system stability, i.e. Vou < 0 the following update laws are proposed: 


ng Aw, Ma; if Gu < GuM 


gu = 4 nz Aw, Mu; if (Gu = guar) &(Aw,M, < 0) (24) 
0; if (Gu = guM)&(Aw,M,, > 0) 


Wy =P Away (25) 
The update rule can be interpreted as: 


= If gu < Gum, it means Gu is within the allowed limit, then ĝu can take any value and here we choose 
Gu = 1, Awu Mau, then Ju (Awu Mu — Nuĝu) = 0. 
— If ĝu = g it means ĝu reaches the upper bound, then Ôu cannot choose a positive value. So if Aw, Mu < 
0 then we will choose ĝu = n7 "Awu Mau < 0, and Ju (Awu Mu — Mugu) = 0. 
— If ĝu = gum and Aw, Mau > 0, we choose ĝu, = 0. This keeps the value g,, equal to the upper bound 
value. Then, ĝu = Ju — Gu = Ju — Jum <0, deduce: Ju (Awu Mu — Nudu) = Ju Awu Mau < 0. 
Thus, in all cases we always have (26): 


Gu(AwuMu a Nugu) < 0 (26) 
In addition, from the update rule (25), we have (27): 


WE (Awuha —T'z!W,) =0 (27) 
Thus: 


Vou = —ku, (Awu)? + WE (Awyhy —Cz!Wa) + Gu(AwuMy — tudu) < -ku (Awu)? <0 (28) 


— Step 3: defining rewind velocity tracking error variables as (29): 


Aw, = Wp — Wr (29) 


Where w, is the reference velocity for the rewind roll. Choose the Lyapunov candidate function: 


V, = (Awr)? = 5 (wr — wra)”. Taking the derivative of V,.: 


V, -i Aw, Aw, = Aw, (Ww, = Wrd) = Aw, (fr + g9rMr z Wrd) (30) 


In order to obtain V, < 0 the control signal M, is chosen as (31h: 
M, _ — 5 (fe — Wra + ky, Awu) (31) 
where kw, is a positive real number, then: 


Vp = —ku, (Awr)? < 0 (32) 


However f, and g, are difficult to determine precisely. Thus, we will approximate them by f and g,, respec- 
tively. Then, the control signal is calculated as 33): 


M, = -4 (fr — Gra + ku, Awr) (33) 


The derivative of V, becomes (34): 
V, = Aw, (fr — i) + fo F (gr a Gr) M, + GrM, — Ùra) 
= Aw, (fr — fr) + Awr (gr — ĝr)Mr, — kw, (Awr)? 


Next, the approximation rules for f, and g, are proposed. Supposing that the function f, can be computed by 
an ideal RBF as (35): 


(34) 
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fa = WEh, (35) 


JT 


|» nis the number of 
T 


neural in the network, and h, is a Gaussian function vector h, = [hp, hro,- pn |° . The approximation 
function f, are calculated using RBF as (36): ; 


Where W, is the ideal weight vector of the neural network W, = [Wri Wio,- We 
A 


f= Wile (36) 


Where W,. = [Wri, Wr, <; Weal? is weight vector. We have to establish the weight updating rule so that the 
system is stable. In addition, we also have to find the approximation rule for gp. It is noted that gr = 1/Jr, 
then there exists lower bound g, > grm > 0. Supposing that g, can be po all by g,. We define 
W, = W, — W,; ; Jr = Gr — Gr. Choose the Lyapunov candidate function V2, as 

Var = = Vr T3 IWIT; W, T ing? = a 3(Aw,)? F IWIT W, T smi? (37) 


where I`, is a positive definite diagonal matrix and n, is a positive real number, then: 


Var = Vp + WITS W, + nõ 
= Awu( fe — fr) + Aws (gr — Gr) My — ku, (Awe)? + WETS W, + ne Grdr 
= Aw WT hy + Awd My — ku, (Awr)? + WIT W, + eGrdr 
= —ku, (Awr)? + WE (Awphy + Po! W,) + Gp(AwuMy + negr) 
= hy, (Aw)? + WE (Awe, — T7 W,) + r (Awr Mr — nedr) 


In order to guarantee system stability, i.e. Vər < 0 the following update laws are proposed: 


Aw, M, rs if g Gr > GrM 
or = Ley ri if (Gr = gra )&(Aw,M, > 0) (38) 
; if (ĝr = gr M)&(Aw,M,. < 0) 


W, =T Awd. (39) 


The update rule (38) can be interpreted as: 

— If g Gr > grm, it means g, is within the allowed limit, then Ôr can take any value and here we choose 
Îr = N7 Aw,M,, then ĝp (Aw, M, — Nrår) = 0. 

— If Gr = grm, it means g, reaches the lower bound, then Gr cannot choose a negative value. So if 
Aw, My > 0 then we will choose g, = n7! Aw, My, > 0, and g,(Aw,M, — mgr) = 0. 

— If ĝ = grm and Aw,M,. < 0, we choose Gr = 0. This keeps the value g,. equal to the lower bound 
value. Then, 9, = gr — Gr = Gr — grm => 0, deduce: g,(Aw,M, — NrGr) = g,Aw,M, < 0. Thus, in 
all cases we always have (40): 

Jr (Awr M, — Gr) < 0 (40) 


In addition, from the update rule (39), we have (41): 


WZ (Awyhp —T=1W,) = 0 (41) 
Thus: 


Vor = —ku, (Awr)? + WE (Aw,h, + 71W) + (Aw M, — må) < ku, (Aw,)? < 0 (42) 


The control laws, as well as the updating rules, have been developed successfully. 
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4. SIMULATION RESULT 


Simulations have been conducted using MATLAB/Simulink to verify the effectiveness of the proposed 
control method. A comparison is also made between the proposed controller and conventional backstepping. 
The simulation parameters are: Ruo = 0.1 m; Rro = 0.05 m; bry = 0.00002533 Nms; br, = 0.00002533 
Nims? Juo = 1.5 kg.m?; Jro = 0.5 kg.m?; a = 0.00005 m; kwu = 230; kwr = 230; and k; = 500. 

In the simulation, for both the conventional BC and the RBF-BC, the reference tension is 10 N. The 
rewind reference speed profile is as follows. The speed increases linearly from 0 to 2 m/s from 0 to 5 second. 
The speed is kept constant from 5 to 20 seconds. Finally, the speed decreases linearly from 20 to 25 second. 
Two situations are considered in the simulation including the system parameters are known precisely, i.e. the 


system parameter error is 0% and there are system parameter error 20%. In both considered situations, the 
simulation results are shown in Figure{a)-(f). 
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Figure 3. Simulation results with no parameter error and 20% parameter error: (a) tension-no parameter error, 
(b) tension-20% parameter error, (c) unwind roll speed-no parameter error, (d) unwind roll speed-20% 
parameter error, (e) rewind roll speed-no parameter error, and (f) rewind roll speed -20% parameter error 


For the case of no error, at times 0, 5, 20, and 25 seconds, due to a sudden change in velocity, there 
are tension overshoots 0.8% as shown in Figure[3{a). After that, the tension returns back to the reference value. 
For the BC controller, the return time is very fast. But for RBF-BC, it takes about 1.2 s to return to the set value 
since the system parameters need to be estimated. In addition, at the steady state, for the RBF-BC, there exists 
a little steady-state error (0.006%). The unwind velocity in case BC and RBF-BC both track the reference 
value; however, in the case of RBF-BC there exists a little error (0.0025%) as shown in Figure Bo). For the 
rewind roll, the velocity tracking is also obtained for both controllers as shown in Figure Bie). When there are 
parameter errors for the BC controller, there is exists a difference between the actual tension and the reference 
value, the difference increases when the error increases. 


For the case of 20% error, the difference is 0.45% as shown in Figure[3{b). For the RBF-BC controller, 
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the difference between the actual and the reference value is much smaller. It is because the BC controller 
depends on the system parameter, but RBF-BC can adapt to parameter variation. This is the advantage of the 
proposed RBF-based controller. For both BC and RBF-BC controllers, velocity tracking is always guaranteed 
as shown in Figures Bid) and (f). This is the advantage of a backstepping-based controller for WTS. 


5. CONCLUSION 

A BC coupled with a RBF neural network is proposed for WTS in this paper. The BC strongly depends 
on system parameters. Thus, the RBF is proposed to update the changing of the system parameter. The weight 
update rules of the neural network have been developed to stabilize the system and guarantee good performance. 
Simulation results demonstrate the effectiveness of the proposed controller compared to the conventional BC. 
In the future, we will continue to improve the performance of the controller to overcome the ”explosion of 
terms” phenomenon in the BC. We also consider implementing the proposed controller in a real system. 
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